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Abstract 

We present an evaluation of the 1-loop prefactor in the lifetime of a metastable 
state which decays at finite temperature by bubble nucleation. Such a state 
is considered in one-component yj 4 -model in three space dimensions. The 
calculation serves as a prototype application of a fast numerical method for 
evaluating the functional determinants that appear in semiclassical approxi- 
mations. 
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1 Introduction 

The decay of metastable states by bubble nucleation appears in a large variety 
of physical contexts. It has received the attention of particle physicists and 
cosmologists due to its possible role in the evolution of the universe. Since 
bubble formation is a basic mechanism in the kinetics of a first-order phase 
transition a precise determination of its rate is of prime importance. The 
semiclassical approach to bubble formation has been developed by Langer 
[|I|, 0] and Coleman and Callan |3|, fj|. The leading factor in the transition 
rate is determined by the classical euclidean trajectory. Quantum corrections 
may however modify the rate in an significant way. Their evaluation for a 
realistic model in three space dimensions is an enterprise that can easily 
reach the limits of practical computability. It is therefore very useful to have 
method that leads to a fast numerical algorithm. 

We will develop here such a method using as a simple model the four- 
dimensional <y9 4 -theory at finite temperature T given by the Euclidean action 



The field potential U((p) is assumed to have two non-degenerate minima 
= and y?+ > (Fig. 1). 

Any state built on the local minimum is metastable. Its decay rate 
per unit volume 7 = T/V at sufficiently high temperature is dominated 
by the energy E = ST of a fluctuation which looks like a bubble <fi{x) of 
the + -phase. This bubble is in unstable equilibrium between collapse and 
unbounded expansion. The tree level approximation determines the order of 
magnitude of the decay rate as 7 m exp{—E((j)(x))/T}. 

Fluctuations around the critical bubble contribute a pre-exponential fac- 
tor to the decay rate which is known to take the form Jl|, [2| 




(1.1) 




(1.2) 



to 1-loop accuracy. The coefficient T> here reads 



V(T) 



det\-(d/dr) 2 - A + U" ((/>)) 
det(-(d/dr) 2 - A + U"{0)) 



(1.3) 
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The prime in the determinant implies omitting of three zero modes in it. 
The temperature dependence of T> arises from imposing periodic boundary 
conditions in the time direction with a period 1/T. T*(T) as introduced in 
Eq. ( |1.3| ) is ill-defined because of ultraviolet divergences. As discussed in [|J 
they are absorbed by the counterterm action S ct which has been introduced 
in the exponent. It will be specified below. 

As we have mentioned above we will present here a fast method for eval- 



uating the fluctuation determinant (1.3). The method which we are going to 



use is based on a well known theorem H that is formulated and proven in an 
elegant way in Appendix A of S. Coleman's 1977 Erice Lectures It was 
also used some time ago in analytical calculations of the determinant (|1.3j ) in 
(l+l)-dimensional space in the thin- wall approximation [|7j, [5| and in Ref. [[]]. 
What is new here is application of that idea to a system in three dimensions, 
the elaboration of a numerical method and the discussion of regularization 
and renormalization. 

We decompose fluctuations around 4>(x) into partial waves, calculate the 
ratio of determinants J\ of radial operators, using the theorem mentioned 
above, and, finally, obtain lnT> as + 1) hi J;. 

In calculating In 22 we exclude the divergent perturbative contributions 
of first and second order in the external field created by the critical bubble. 
The regularized values of these contributions are then added analytically. All 
divergences of In T> appear in the standard zero-temperature tadpole and fish 
diagrams. 

This paper is organized as follows: In the next section we specify the 
form of the potential, write the equation for the critical bubble and present 
our numerical results for E(<p). In Sec. |3] we describe the calculation of the 
regularized fluctuation determinant (|1.3| ). A possible renormalization scheme 
is applied to the result in the Sec. f§. 



2 The Tree-Level Energy 

In this section we discuss classical properties of the critical bubble. The 
generic one-component 4 -potential reads 

[%) = imV~W 3 + ^ 4 (2.1) 
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We choose the same dimensionless variables as in Ref. [|IU| , namely x = X/m, 
t = u/m, <p = The energy of a time-independent fluctuation then takes 

the form 

m = %\^ {I (V$) 2 + - \^ + ^ 4 ) • (2.2) 

where 

Am 2 

a = ^. (2.3) 

The limit a — > 1 corresponds to the thin wall approximation. 

The critical bubble is a spherically symmetrical stationary point of E 
( |2.2|) obeying a dimensionless Euler - Lagrange equation 

d 2 $ 2rf$ _ 3 o 1 3 n , 0/1 . 

— - H $ + -$ 2 $ = . (2.4) 

dR 2 RdR 2 2 V ; 

We have solved this equation using the shooting method. The profiles $(-R) 
are shown in Fig. 2 for some values of the parameter a . 

We parametrize the value of E in the same way as in Ref. flOU : 

4.851 m 3 „, , . 

E = f(a) 2.5 

rj 2 

The function /(a) is plotted in Fig. 3. 



3 Calculation of the Fluctuation Determinant 

In this section we discuss a method of computing the ratio of functional 
determinants (|1.3| ). which is based on earlier papers 0-|J. 
The explicit form of the operator in the nominator ( |1.3| ) is 

-7T-2 - A + m 2 + V(r) (3.1) 

with the periodical boundary conditions for the eigenf unctions \I/(t, x) = 
*f?(r+ 1/T,x). The space-dependent part V Q3.1J ) of U"((f)) (|1.3|) reads 

7( r ) = [/"($) - m 2 = -67$(r) + ^A0 2 (r) . (3.2) 
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The free operator in the denominator ( |1.3| ) takes the same form as ( |3.1| ), but 
with V{r) = 0. 

The time-independence and spherical symmetry of the background bubble 
field yield a classification of uj^ 1t1t , the eigenvalues of (|3TT|), with respect to 
the number n of their Matsubara frequencies v = 2irnT, radial quantum 
number n r , and angular momentum I. One can formally write the ratio of 
determinants (|1.3| ) as 



oo oo oo 



v<?) = n n n 

n=— oo 1=0 n r =0 



^n,l,n r 
n,l,n r 



21+1 



(3.3) 



M") = —T^y (3-4) 



with uj^]„ standing for the free eigenvalues. 

The outline of our calculation is as follows. We compute first (Sec. |3.1|) 
for each partial wave the product over n r , i.e. the ratio of the determinants 
of the radial operators 

detH v i 
detH® 
here 

H vJL = --f- 2 ---^ + " 2 + ™ 2 + V(r) (3.5) 
ar A r ar 

and lif} has the same form, but without V(r). 
As the next step we calculate the product over /: 

oo 

F(i/)=£(2J + l)lnJ,(i/) (3.6) 

1=0 

(section |3.2|) The function F(u) is the sum of all three-dimensional one-loop 
one-particle-irreducible diagrams. 

In terms of this function the final result reads 

oo 

lnP(T) = Y, F(2ixnT) . (3.7) 

n=— oo 

This expression is formal because of the ultraviolet divergences in it. We first 
evaluate ln£> without the divergent part which is then accounted for in Sec. 
IO. 
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3.1 Determinants of the Radial Operators 

In order to find Ji(v) ( |3.4| ) we make use of a known theorem |5|, |6| whose 
statement is 

detH Ujl ^v,i{r) , . 
j^r = hm — ^ . (3.8) 

detHf) ™VS } W 
Here ip v j and tp^j are solutions to equations 

H U ^ V>1 = , <M° } = (3.9) 

and have same regular behavior at r = 0. More exactly, the boundary 
conditions at r = must be chosen in such a way that the right-hand side of 
Eq. (|3.8|) tends to 1 at v — > oo. 

It is convenient (Tl| to introduce a function h{r) such as 



when ip^i is chosen here to be a spherical Bessel function 



^,i = (l + /ii(r))i t (/cr) , fc(0) = (3.10) 
3 to be 

i,(«r) = f— V /2 /,, x(«r) , k(z/) 2 = z/ 2 + m 2 . (3.11) 

Therefore, by the theorem ( |3.8j ), the ratio of determinants ( |3.4j ) can be ex- 
pressed as 

J t {u) = (1 + hi{oo)) . (3.12) 
In terms of the h function the first equation (|3.9|) reads 



(-^ + 2(^ m + 1 -)±)h l (r) = V(r)(l + h l (r)) . (3.13) 
y ar z \ti{Kr) r J ar J 

It is worth to consider the structure of a perturbation expansion 

oo 

h{r) = Y.h?\r) (3.14) 

k=l 

in powers of the potential V(r). This entails an analogous expansion for the 
ratios Ji(v) in the sense that j\ k ^ = h\ k \oo). The fc-order contribution h\ k ^ 
obeys an equation 

Lf 2+2 (^ m + \)f\ h f\ r)=V(r) ^ (r) , (3.15) 
V ar z \ti{Kr) r I ar I 
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hf ] = 1. The same equation is valid when h\ are replaced by h\ 

J2^Lk hi ■ In this notation hi = 
to equation (|3.15|) in the form 



Yl,T=k h\ q ^ ■ In this notation hi = h\ . A Green function that gives the solution 



h9'\r) = - dr'r' 2 G l (ry)V(r f )h?~ 1) {7 J ) (3.16) 
Jo 

with the right boundary condition at r = reads 

G^r') = k {iii^hiKT^^- -kiKr'MKr')] . (3.17) 

Here r< = min{r, r'}, r> = max{r, r'} and 

/ 2 \ 1/2 

*;,(*) = — Ki(z) . (3.18) 

The first term on the right-hand side of Eq. (|3.17 ) does not contribute to 
hf : \oo). The Green function (|3.17| ) gives rise to connected graphs as well as 
disconnected ones (Fig. 4). The latter are cancelled in ln(l + /^(oo)) whose 
expansion in fc-order connected graphs ./zcon(ft) reads 

oo /_i \k+l 

In J t (y) = ln(l + h^oc)) = £ J*?Uv) . (3.19) 

k=i K 

This formula is analogous to the expansion of the full functional determi- 
nant in terms of Feynman diagrams 



ln£>(T) = ^^-f A {k \T) . (3.20) 



Here A^(T) is the 1-loop Feynman graph of order k in the external potential 
V(\x\). 

Indeed, it is obvious from Eq.( |3.16| ) that h\ and, therefore, j\ k l on are 
of the order V k . Since the expansion of lnX> in powers of V is unique, we 
conclude that 



A«(T)= £ £(2/+l)j/tL(2™T) . (3.21) 

n=— oo Z=0 
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One can verify this relation explicitly by expanding the propagator in 
as 



/ 



d 3 p 



(2vr) 3 p 2 



m 2 + v 1 



oo 

K£(2Z + l)i,(K|f|)A;,(«|jfl)fi 

1=0 



x ■ y 

W\\y\. 



(3.22) 



and performing the integration over all angular variables in the x-represen- 
t at ion. 

It is not difficult to check - making use of a uniform asymptotic expansion 
of the modified Bessel functions in (|3.17| ) - that jw \j\ k as / — » oo. That 
results in the expected quadratic and logarithmic ultraviolet divergences in 
InD due to the contribution of h^(oo) and h^ 2 '(oo). We have computed 
numerically lnD® which is the sum ( |3.20| ) without first and second order 
diagrams and A^ 2 \ It reads explicitly 



2>W(T)= ^ F^(2nnT)= £ £(2Z + 1) (In Jfr)) 



(3) 



(3.23) 



-oo 1=0 



where 



ln(l + fy(oo)) -^ (1) (oo) 



fcf 2) (oo) 



r /I 



OO 



(3.24) 

(2) 

The terms in square brackets here correspond to the fish diagram J\ con (Fig. 
4). Since all contributions to lnD^ are ultraviolet finite, we need no regu- 
larization in computing them. The divergent contributions of the first and 

second order in V will be considered in Sec. Ol. 

We have determined hi(r) as solutions of Eq. fl3.13|) , and h\ (r) as that of 



hf\oo) have 



Eq. (|3.15|) by Nystrom method. The values of h\ (oo), h\ (oo 
been evaluated by performing integration ( |3.16| ). Only the last term in the 

Green function fl3.17| ) contributes here since r — > oo. The ratio of hf^^oo) 
found via differential equation ( 3.1 5|) to that calculated as the integral ( [3.16| ) 
has been used to control the accuracy. In order to avoid numerical subtraction 
that might be delicate we re-write the term ( |3.24| ) to be summed up on the 
right-hand side ( |3.23p in the form 



(lnJ^)) (3) 



ln(l + fy(oo)) - h{oo) + -h^oo) 2 



+ h 



(3), 



1 



oo 



hf\oo) 



hAoo) + h, (oo 



(3.25) 
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Each of the three terms on the r.h.s. is now manifestly of order V 3 . The 
subtraction done in the square bracket is exact enough when the logarithm 
is calculated with double precision. 

In the numerical computation /i/(oo) is to be replaced, of course, by 
hi{r max ). We have found that hi(r) becomes constant within relative de- 
viation of O(10 -6 ) for r « (12 18)/m and we have chosen r max in this 
range of values. 

We have neglected till now the existence of the negative mode < 
and three zero modes u;q 10 = 0. The former results in negative value of 
Jq{v) = 1 + h (oo) at v = 0. According to Eq.(|1.2j) one has to replace 



to, 



0,0,0 



by |co>o 00 |. That implies taking the absolute value of Jo(0) in Eq. ( |3.23| ). 

The zero modes manifest themselves by the vanishing of Ji(z^) at v = 0. 
In the numerical calculation this zero was found at Uq x = O(10 _5 m 2 ). 

It can be easily seen |7], |8|, |J that exclusion of the zero modes implies 
replacing of Ji(0) by its derivative 



d 

dp 



J'M = -JTTT^hi (oo)U . (3.26) 



In Fig. 5 we present some results for the functions hi. The values of the 
first h^\oo) are plotted vs. 21 + 1. For the terms summed in Eq.( p.23[) we 
have found good agreement with the expected behavior 1/(21 + l) 4 . 

3.2 Calculation of 



Our next step is performing summation over / in Eq. (|3.23| ). It has been done 



by cutting the sum at some value l max and adding the rest sum from l max + 1 
to oo of terms fitted with 

, Const Const' 

lnJ < ~WTW + WTW ■ (3 ' 27) 

The summation was stopped when increasing of l max by unity did not change 
the result within some given accuracy 5. This happened, for example, at 
lmax — 12 at v = 0, a = 0.5. The value of 5 was automatically increased 
during the calculation if the corresponding accuracy had not been reached. 

The convergence becomes worse at higher v or a. The reason is that 
the asymptotic behavior ( |3.27| ) sets in at / ^> [y 2 + m 2 ) 1//2 r e // where r e ff is 
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the characteristic size of the bubble. It is of order 1/m at small values of a 
and can be estimated as (4/ (3(1 — a)m) + 0(l))m _1 near the thin-wall limit 
a — > 1. As the maximal value of the angular momentum that we have used 
is / = 30, the computations have been stopped somewhere at v ~ 10m. 

The resulting F^'{y) is shown in Figs. 6-8. Its magnitude at v — 
gives the value of the infinite-temperature determinant ratio lnP^oo). To 
illustrate the efficiency of the method we note that the evaluation of F^{v) 
for one value of v takes typically 10 -j- 30 sec CPU time on a standard PC 
with 486 processor. 

The finite temperature is accounted for in lnP^ 3 ^ which is computed ac- 
cordingly to Eq.(|3)- A s F®(v) ~ 1 /z/ 3 at high z/, the summation over 
Matsubara frequencies is elementary. 



3.3 Inclusion of and 

We have found the value ln2)( 3 ) which is the sum of all 1-loop diagrams of 
the third order and higher. Now we add to the result first- and second-order 
finite-temperature Feynman graphs (T) and (T) ( |3.20| ) calculated in 
the standard technique. It is convenient to represent both of them if the form 

A (fc) (T) = [A {k) (T) - A {k) (0)] + A (fc) (0) (3.28) 

where the values in square brackets follows. 

A W(T) - A«(0) = i jT drr 2 V(r)Q (3.29) 

with 

r°° / \ 1/2 
Q{z) = dy{l + y 2 ) [coth(zy) - 1] . (3.30) 

The second order gives an UV-finite contribution to F(u): 

( 2 \ "I/ 2 2 
1+f — j j (^J drrV (r) sin qr^j 

(3.31) 

We have taken A (2) (T) - A (2) (0) ( [3.28D into account by calculating nu- 
merically the difference 

A^ 2 \T)-A^(0)= V F^(2nnT)-- f° — F®{v) . (3.32) 

„_ „ T J-oo 27T 
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It is sufficient to make an intermediate cutoff regularization in order to de- 
termine uniquely the difference betwen these two logarithmically divergent 
quantities. 

All UV divergences have been moved now to the last terms in (|3.28 ) which 
are the standard tadpole and fish diagrams at zero temperature. 

To sum up, we have calculated the functional determinant ( |1.3| ) as a sum 
of following contributions. 

In (m 6 P(T)) = [A {1) (T) - A {1) (0)] + [A (2) (T) - A {2) (0)" 

oo 

+ F^(27inT) + A^(0) + A^(0) . (3.33) 



Here the first term is given by Eq.( [3.29| ), the second one is defined in ( p. 32 



3.31| ), the sum of has been calculated numerically, and two last terms are 
usual zero-temperature Feynman diagrams ( 3. 20]) in the external potential 
V. 

AW(0) and A^(0) contain now all the ultraviolet divergences and have 
to be regularized. The cutoff dependence introduced thereby disappears, 
however, in the full 1-loop contribution to the effective action 

Si-ioo P = \ Hm 6 V(T)) + S ct (3.34) 



that enters the formula for the transition rate Ql.2| ). In the model considered 
here S ct has the form 

S ct = J d 4 x Ltp + -5m 2 <p 2 - 5riv 3 + ^ V) (3-35) 

A possible scheme for fixing the counterterms is deferred to the next section. 



4 A Possible Renormalization 

While renormalization requires just standard techniques it is not straightfor- 
ward here to select a specific renormalization precription because the scheme 
depends strongly on the physical context in which the first-order phase tran- 
sition is considered. Though we consider Eq. Q3.33j ) in its general form as 
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our main result, we would like to discuss now a possible scheme of fixing the 
counterterms. 

The scheme is chosen in spirit of renormalization suitable in consideration 
of the electroweak cosmological phase transition. As the latter is tightly 
connected with particle physics, it is appropriate to fix the mass and the 
vacuum expectation value of the scalar field. This implies the following 
conditions on the temperature-dependent effective potential U e ff((p;T) and 
the 1-loop Euclidean propagator in the true vacuum D + (p 2 ): 

DA-mV) = 



U' eff (<l> + ;0) = 

U eff(<f>+ 



0) = 

0) = U(<P + ) (4.1) 



where m 2 + stands for the particle mass in the true vacuum. The last two 
conditions set the false vacuum to be at (p — and fix the density of energy 
stored in it. They are more specific for our toy model. 

In addition to Eqs. (|4.1| ) one has to solve one more system of equations 



I m V+-#+ + ~A0i = U{<f>+) 
m 2 - 3r)<f> + + ~\(f>l = 

m 2 - 6r]<f) + + -X(f)l = (4.2) 

in order to express m 2 , rj and A in the classical bubble energy ( |2.2j ) ,( [2.3| ) in 
terms of m+, <p+ and U((f)+). 

The results of the application of this renormalization scheme are plotted in 
Figs. 9-11. We note that the temperature appears neither in the field poten- 
tial ( [2.1|) nor in the renormalization conditions ( |4.1|) . All finite-temperature 
corrections are therefore contained in the fluctuation determinant. This re- 
sults in a linear temperature dependence due to the first term in (|3.33|) . 



Thus, our simple 1-loop approximation fails at very high temperature. It is 
not valid also at a — > due to the high difference of the mass scales of states 
built on the false and true vacua. This manifests itself by a logarithmically 
large contribution appearing in (|3.35|). 
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5 Discussion and Conclusion 



The model considered here is only semi-realistic. Nevertheless we would like 
to add some remarks on our results. We find that the correction to the 
bubble nucleation rate ( |1.2|) coming from diagrams of the third order and 
higher favors the transition (i.e. InT)^ < 0). The sign of the full 1-loop 
contribution to the effective action Si-i oop depends on the renormalization 
scheme. With our choice ( f4.1|) it enhances j(T) at high temperature due to 
the tadpole diagram ( |3.29| ). This contribution becomes too big at T ^> m and 
then the naive application of the 1-loop approximation becomes inconsistent. 
This is a manifestation of the known problem of relating the parameters of 
the theory at zero and very high temperatures. Another feature of our results 
(Fig. 9-11) is a weak logarithmic singularity at T = T min , the temperature 
at which one more fluctuation mode of the critical bubble becomes unstable. 
Numerically it is unimportant. Moreover, in a finite region of temperature 
Tmin < T < T tunn quantum tunneling has to be taken into consideration. 

We have developed here a method for calculating functional determinants 
and have tested it under realistic conditions. In particular we have shown 
that within this formalism the problems connected with ultraviolet diver- 
gences, zero and unstable modes can be handled easily and without loss of 
numerical accuracy. Moreover the algorithm is so fast that the whole calcula- 
tion presented here can be performed on a standard PC with a 486 processor 
within one hour for one value of a. The extension of the method to more 
complex gauge and fermionic systems is straightforward once the fluctuation 
equations are known [O, [13]. 
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Figure Captions 



Fig. 1 Potential U(<p) fl2~T| ). It is plotted in dimensionless form which enters 
the integral ( |2.2|) . The curves are labeled with the value of a. 



Fig. 2 The bubble-profile functions <&(R) in units defined after Eq. ([2.1|) at 
a = 0.1, 0.5, 0.8, 0.9 . The radius of the thin- walled bubble mr tw ~ 
3/(4(1 — a)) is marked with a dashed line for a = 0.9 . 

Fig. 3 The values of f(a) (j2J^) (solid line) and -lOu; 2 /m 2 (dotted line) vs. 
a. 

Fig. 4 The structure of first h\ (oo). The solid line represents the last term 
in the Green function (|3.17| ). Dots stand for V(r). 



Fig. 5 The values of j/^, n at A; = 1 (curve 1), k — 2 (curve 2) and jfl on 
(curve 3) against 21 + 1 in double logarithmic scale at a = and v = 0. 
The straight-line behavior at large I corresponds to the expected power 
law (2Z + l)~ 2fc+1 . 

Fig. 6 Absolute value of F^(u) at a = 0.1 vs. 21 + 1 in double logarithmic 
scale. A logarithmic singularity is seen at v = The actual value of 
F^{v) is negative above this point. In the region v < oj- it is complex 
and has no physical meaning. The dotted line represents expected 
relative error in F^ 3 \u) which is estimated as 56 (see subsection |3.2| ). 

Fig. 7 The same as in Fig. 6 at a = 0.5 . 

Fig. 8 The same as in Fig. 6 at a = 0.8 . 

Fig. 9 The value of Si^i oop ( |3.34| ) at a — 0.1 vs. T/m: The curves 1-3 
correspond to the first, second and third terms in r.h.s. (|3733| ) . The 



fourth one represents the sum of last two terms in ( |3.33| ) and S ct ( |3.35| , 
|4.1|) . This contribution depends on temperature via the factor 1/T 
only. Curve 5 displays the full result Si-i oop . The temperature range 
is bounded from the left by T min = 0.2301m where the result has a 
logarithmic singularity. 

Fig. 10 The same as in Fig. 9 at a = 0.5 ; here T min = 0.148m. 
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Fig. 11 The same as in Fig. 9 at a = 0.8 ; here T min = 0.0598m. 
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